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Two  numerical  approaches  to  the  problem  of  theoreTfcal ly  predicting 
closure  pulses  produced  by  bulk  cavitation  are  reviewed:  a  single-stage, 
finite-difference  approach  based  directly  on  the  governing  differential 
equations  of  hydrodynamics  and  a  multi-stage  approach  which  assumes  several 
different  simplified  formulations  of  the  governing  equations  at  different 
times  during  the  bulk  cavitation  process.  New  extensions  to  both  current 
approaches  are  reported  by  the  author  and  used  to  solve  an  example  problem. 

In  particular,  a  new  graphical  method  for  predicting  the  form  of  c  osure 
pulses  everywhere  except  near  regions  of  subsonic  closure  is  described. 

The  two  numerical  approaches  are  compared  on  the  basis  of  the  example 
problem  solutions  and  several  aspects  of  the  methods  which  require  further 
development  are  identified.  l~j  *| 
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NUMERICAL  PREDICTION  OF  SUPERSONIC  BULK-CAVITATION  CLOSURE  PULSES 

by 

Joseph  A.  Clark 

Acousto-Optics  Laboratory,  Catholic  University  of  America 
Washington,  D.C. ,  20064,  U.S.A. 

I.  INTRODUCTION 

Bulk  cavitation  is  produced  by  the  interaction  of  shock  waves  from 
underwater  explosions  with  the  water  surface.  As  a  surface  layer  of  water 
is  accelerated  upwards  by  the  reflecting  shock  waves,  a  region  of  bubble- 
filled  water  forms  below  it.  When  the  surface  layer  falls  back,  the 
cavitation  region  rapidly  collapses  and  a  second  pressure  pulse  is 
produced  as  it  finally  closes.  This  closure  pulse  associated  with 
bulk  cavitation  is  of  practical  significance  because  the  energy  flux  density 
(time  integral  of  instantaneous  intensity)  of  the  closure  pulse  is  often 
comparable  to  that  of  the  initial  shock  wave. 

Characteristic  features  of  bulk  cavitation  closure  pulses  can 
be  noted  in  the  example  shown  in  Figure  1.  This  example  shows  the  response 
of  a  pressure  gage  located  70  feet  below  the  surface  and  500  feet  away 
from  a  charge  of  1200  lbs  of  HBX-1  explosive,  which  was  detonated  70  feet 
below  the  surface.  (The  pressure  history  was  taken  from  data  collected  by 
the  Underwater  Explosions  Research  Division  of  David  Taylor  Naval  Ship 
Research  and  Development  Center.)  The  initial  shock  wave  (A)  passes  by 
the  pressure  gage  until  it  is  cutoff  by  a  surface  reflected  shock  wave  at 
(B) .  Then  the  gage  pressure  drops  to  a  constant  negative  value  (corresponding 
to  zero  absolute  pressure).  This  condition  indicates  that  the  gage  is  near 
a  region  of  cavitating  water  which  collapses  at  about  (C) .  A  closure  pulse 
then  passes  the  gage  which  is  of  lower  peak  amplitude  but  much  longer  duration 
than  that  of  the  initial  shock  wave  pulse.  Other  examples  of  closure  pulses 
have  appeared  in  the  literature  [1-3]. 

Theoretical  methods  for  predicting  bulk  cavitation  closure  pulse 
characteristics  have  followed  two  basic  approaches  in  their  development.  Either 
a  single  hydrodynamic  formulation  of  the  entire  bulk  cavitation  problem  has 
been  attempted  or  else  various  stages  of  the  problem  (propagation  of  the 
initial  shock  wave,  formation  and  collapse  of  the  cavitation  region,  prop- 


(4) 


agation  of  the  closure  pulse,  reformation  of  the  cavitation  region 
have  been  formulated  separately.  This  investigation  originated  in  an  attempt 
to  compare  theoretical  predictions  of  bulk  cavitation  phenomena  with  experi¬ 
mental  results  obtained  by  producing  bulk  cavitation  in  a  laboratory 
facility  [4,5].  It  was  found  that  several  extensions  of  the  available 
theoretical  prediction  methods  were  needed  before  useful  comparisons  between 
theory  and  experiment  could  be  made. 

In  this  paper  the  present  state  of  development  of  theoretical  bulk 
cavitation  prediction  methods  will  be  assessed.  Prior  work  will  be  briefly 
reviewed.  Recent  extensions  of  the  methods  will  be  described.  Theoretical 
results  obtained  for  a  typical  problem  will  be  presented  as  a  basis  for 
comparing  the  two  approaches.  Also  several  aspects  of  the  methods  which  require 
further  development  will  be  identified. 
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2.  HYDRODYNAMIC  APPROACHES 

Bleich  and  Sandler,  in  one  hydrodynamic  approach,  numerically  solved 
a  one-dimensional  bulk  cavitation  problem  by  the  method  of  character i sties [6] . 
They  assumed  linearized  equations  of  motion  and  continuity,  together  with 
a  bilinear  equation  of  state  [7-93  *  Sandler  later  developed  a  finite-difference 
method  for  solving  the  same  formulation  of  the  1-D  problem  which  gave  results 
directly  as  pressure,  velocity  or  density  histories  that  could  be  compared 
with  experimental  measurements.  More  recently  he  has  developed  a  two- 
dimensional,  finite-difference  method  [10], 

Following  suggestions  by  Sandler,  this  author  extended  the  1-D,  finite- 
difference  method  to  allow  for  the  specification  of  arbitrarily  shaped  initial 
shock  waves  and  a  variety  of  bottom  boundary  conditions  (rigid,  soft,  time- 
varying  ...).  Results  obtained  by  a  Fortran  impl imentat ion  of  this  method  will 
be  described  later  in  this  paper. 

3.  MULTI-STAGE  APPROACHES 

Step-wise  numerical  integration  of  the  governing  differential  equations 
can  be  avoided  until  cavitation  occurs  If  an  analytical  model  of  the  explosion 
shock  wave  proposed  by  Weston  [II]. is  assumed: 

P(r,t)  -  A*((W1/3)/r) -exp(-t/e) 

where  A,B  and  0  are  empirically  determined  constants,  W  is  the  charge  weight, 
r  is  the  distance  from  the  explosion  source  to  the  point  of  interest  and  t  Is 
the  time  after  the  shock  front  arrives  at  the  point  of  interest. 
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Arons  et  al  [12]  showed  that  the  boundaries  of  the  cavitation  region 
could  be  determined  by  three  additional  assumptions.  One  assumption 
expresses  the  bilinear  nature  of  water.  The  second  replaces  the  surface- 
reflected  wave  with  a  negative  pressure  wave  from  an  image  charge  source 
above  the  surface.  The  strength  of  this  image  source  is  assumed  to  become 
weaker  as  the  wave  propagates  in  a  manner  which  just  accounts  for  the 
canceling  of  the  initial  shock  wave  pressure  in  cavitation  regions.  The 
third  assumption  locates  the  bottom  boundary  of  the  cavitation  region 
at  the  depth  where  the  image  source  strength  stops  decreasing.  With  these 
assumptions,  the  upper  and  lower  cavitation  boundaries  at  a  given  horiz¬ 
ontal  range  can  each  be  found  by  solving  a  single  differential  equation. 

Aron's  approach  has  been  employed  by  several  later  investigators  [13- 
15,1]  to  determine  the  extent  of  the  region  of  cavitation  produced  by 
underwater  blasts  of  a  given  charge  weight  and  depth.  The  "kickoff" 
velocity  of  water  particles  within  the  cavitation  zone  just  after  cavita¬ 
tion  begins  can  also  be  determined  by  this  approach.  That  information 
completes  the  solution  of  the  first  stage  of  the  bulk  cavitation  problem. 

Cavitation  boundaries  and  kickoff  velocities  serve  as  initial  condi¬ 
tions  for  the  second  stage  of  the  cavitation  problem.  Costanzo  and  Gordon 
[15]  have  developed  a  numerical,  finite-element  method  for  solving  the 
governing  equations  during  this  stage.  They  simplify  the  equations  to 
exclude  horizontal  motions  and  compressibility  effects.  No  coupling  is 
assumed  to  exist  between  elements  until  an  impact  occurs  with  neighboring 
elements.  Then  rigid  coupling  is  assumed.  The  solution  to  this  stage 
of  the  cavitation  problem  is  given  in  terms  of  the  depth,  time  and  (relative) 
velocity  of  the  last  elements  to  impact  in  the  cavitation  zone  at  a  given 


range. 
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An  Impl imentation  in  Fortran  code  of  the  multi-stage  method  of 
Costanzo  and  Gordon  for  determining  closure  time,  depth  and  velocity 
was  prepared  by  this  author.  This  imp  liman  tat  ion  also  Incorporated  one 
extension  of  the  method  which  accounts  for  the  retardation  of  the 
bottom  boundary  motion  by  impacts  with  elements  of  the  cavitating 
water.  This  extension  was  suggested  by  A.  Misovec  and  developed  in 
close  collaboration  with  Costanzo  and  Gordon. 

Closure  time,  depth  and  velt.-ty  data  provide  the  initial  conditions 
for  the  third  stage  problem  of  predicting  the  form  of  the  closure  pulse. 
Only  the  recent  paper  by  Costanzo  and  Gordon J 5]  has  previously  attempted 
to  predict  the  form  of  the  closure  pulse  by  the  multi-stage  approach  and 
their  solution  only  applies  at  the  single  horizontal  range  where  closure 
first  occurs.  Furthermore,  they  assume  that  the  closure  wave  is  plane 
and  propagates  vertically  in  this  region  of  first  closure.  That  assumption 
is  not  in  general  valid  because  the  closure  path  is  usually  not  horizontal 
in  the  region  of  first  closure.  However,  their  approach  might  give  a 
useful  first  order  approximation  to  the  actual  closure  pulse  form. 

In  order  to  obtain  a  more  realistic  prediction  of  the  closure  pulse 
form,  a  graphical  method  which  accounts  for  the  actual  obliquity  and 
curvature  of  the  closure  wave  was  developed  by  this  investigator.  The 
method  assumes  closure  time,  depth  and  velocity  data  as  initial  conditions. 
It  is  applicable  in  regions  where  the  velocity  of  the  point  of  closure 
along  the  closure  path  achieves  supersonic  speeds.  A  description  of  this 
graphical  method  for  solving  the  third-stage  of  the  bulk  cavitation 
problem  will  be  included  in  the  discussion  of  an  example  problem  in  the 
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next  section. 
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4.  EXAMPLE  PROBLEM 

The  response  in  water  to  40000  lbs.  of  HBX  explosive  detonated  at 
a  depth  of  200  feet  will  be  used  to  compare  results  obtained  by  the  two 
approaches  discussed  above. 

Data  obtained  by  the  1 -D, f in? te-d i f ference  method  can  predict  dynamic 
pressure,  particle  velocity  and  density  histories  in  the  column  of  water 
directly  above  the  explosion.  Figure  2  shows  the  surface  velocity  history 
and  ten  pressure  histories  at  depths  from  10  to  100  feet.  The  propagation 
upwards  of  the  initial  shock  wave  is  observed  in  each  of  the  records. 
(Oscillations  in  the  height  of  the  wave  peak  from  record  to  record  are 
effects  of  sampling  errors  in  the  graphics  program.)  As  the  wave 
approaches  the  surface,  it  is  seen  to  lose  its  exponentially  decaying  tail 
and  instead  drops  almost  instantaneously  to  a  negative  pressure.  This  is 
the  result  of  cancellation  by  the  surface  reflected  wave,  which  occurs 
without  any  additional  assumptions  in  this  approach.  The  initial  shock 
wave  in  the  20-ft.  depth  record  exhibits  a  form  similar  to  that  of  the 
wave  in  the  measured  data  of  Figure  1.  The  surface  velocity  is  seen  to 
increase  and  decrease  rapidly  at  first,  until  cavitation  begins  in  the 
water  beneath  it.  After  that  time,  the  surface  layer  velocity  decreases 
more  slowly.  It  continues  to  decrease  to  negative  values  (which  corresponds 
to  the  surface  layer  falling  back  down  on  the  water  below),  until  the 
closure  pulse  arrives  at  the  surface.  At  that  time  the  surface  velocity 
rises  rapidly  to  zero.  (This  rapid  deceleration  of  the  surface  layer 
would  be  transmitted  to  any  vessel  floating  in  the  water.) 

Closure  occurs  first  near  the  40-ft.  depth.  This  is  indicated  in 
Figure  2  by  the  absence  of  "precursor"  waves  in  front  of  the  closure  pulse 
at  that  depth.  The  precursor  waves  are  associated  with  impacts  between 
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the  cavltating  water  elements  and  the  upper  and  lower  layers  of  water 
which  is  no  longer  cavitating.  The  form  of  the  closure  pulse  Is  approx¬ 
imately  rectangular.  The  duration  of  the  closure  pulse  is  very  short 
near  the  surface.  It  increases  with  depth  until  the  closure  depth  is 
reached.  At  greater  depths,  the  duration  remains  approximately  constant. 

The  predicted  magnitude  and  duration  of  the  closure  pulse  at  the  40-ft 
depth  are  403  ps  i .  and  1**  milliseconds. 

A  new  feature  of  this  finite-difference  program,  as  noted  above,  is 
that  a  variety  of  boundary  conditions  can  also  be  specified.  For  example, 
if  the  original  problem  is  modified  to  include  a  rigid  boundary  at  a 
depth  of  900  ft.  the  predicted  response  of  the  water  is  that  shown  in 
Figure  3.  In  this  case  a  second  shock  wave  propagates  up  into  the 
cavitation  region  about  300  milliseconds  after  the  first.  This  is  the 
wave  reflected  from  the  bottom.  It  does  not  reach  the  water  surface  but 
instead  is  reflected  by  the  cavitating  water.  The  figure  shows  that  the 
second  pressure  pulse  causes  the  closure  pulse  to  be  generated  14  milli¬ 
seconds  earlier.  The  magnitude  of  the  closure  pulse  is  reduced  by  about  to 
385  psi. 

Major  qualitative  features  of  the  bulk  cavitation  closure  pulse  are 
seen  to  be  revealed  clearly  and  directly  in  data  obtained  by  this  finite- 
difference  approach.  For  example,  the  unusual  form  of  a  surface-reflectec 
wave  gradually  attenuating  with  depth,  which  must  be  treated  as  an  ad  hoc 
assumption  in  Aron's  approach,  occurs  automatically  in  this  case.  Also 
complications  of  the  problem  which  often  occur  in  actual  underwater 
b1asts[l],  such  as  reflections  from  the  bottom  or  interactions  with  a 
vessel  floating  on  the  water,  can  be  accounted  for.  However,  these 
solutions  require«very  fine  computat ional  grids.  In  the  present  examples, 
grids  of  120  x  4000  points  were  used  and  the  results  were  still  some- 


what  sensitive  to  grid  size.  (An  increase  in  the  grid  size  to  160x4000 
points  caused  the  predicted  magnitude  of  the  closure  pulse  to  increase  by 
12%  from  403  to  458  psi.)  Also,  the  1-D  code  results  described  above  might 
not  be  directly  applicable  to  at-sea  tests  because  effects  of  pulsing  and 
rising  of  the  explosion  bubble  produced  by  underwater  blasts  have  been 
neglected  [16]. 

The  1-D  finite-difference  method  was  used  in  the  above  example  rather 
than  a  2-D  version  because  comparisons  with  laboratory  experiments  require 
a  computer  code  which  can  be  run  many  times  to  investigate  the  dependence 
of  theoretical  results  on  various  parameters.  The  l-D  code  is  already  just 
marginally  acceptable  in  this  respect.  The  computational  requirements  of 
the  2-D  codt  ire  so  large  that  it  has  been  run  only  a  few  times  on  the 
largest  government  computers  and  results  obtained  with  this  2-D  code  have 
not  reproduced  fine  details  of  the  closure  pulse  form  [17].  The  2-D  code 
in  its  present  state  must  therefore  be  judged  unacceptable  as  a  theoretical 
tool  for  parametric  studies  of  bulk  cavitation  mechanisms. 

Two-dimensional  multi-stage  methods  present  a  very  different,  but 
complementary  view  of  bulk  cavitation  phenomena.  Solutions  of  the  first 
stage  problem  determine  maximum  boundaries  of  the  region  near  an  underwater 
blast  where  water  cavitates  at  any  time  during  the  period  of  interest.  The 
depth  of  the  upper  and  lower  boundaries  of  this  region  are  plotted  as  a 
function  of  horizontal  range  in  Figure  4.  (Data  used  to  plot  the  graph  are 
given  in  Appendix  A.)  The  cavitation  zone  in  this  example  extends  almost 
from  the  surface  to  a  depth  of  50.1  ft.  in  the  region  directly  above  the 
blast.  At  larger  horizontal  ranges  its  depth  increases  until  at  a  range  of 
1050  ft.  its  upper  and  lower  boundaries  are  at  1.8  and  112.6  ft.  respectively. 

m 

The  cavitation  zone  decreases  in  depth  with  further  increases  in  range,  out 
to  a  maximum  range  of  about  2825  ft. 
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Solution  of  the  second  stage  of  the  bulk  cavitation  problem  identifies 
the  path  of  closure  and  final  impact  velocity  as  functions  of  range.  These 
data  are  also  shown  in  Figure  4.  This  solution  predicts  that  closure 
occurs  at  a  depth  of  12.2  ft.  and  at  a  velocity  of  42.3  ft. /sec.  directly 
above  the  blast.  The  predicted  closure  depth  increases  to  a  maximum  depth 
of  68.5  ft.  at  a  horizontal  range  of  1800  ft.  The  impact  velocity  is  0.6 
ft. /sec.  at  that  range. 

The  hydrodynamic  pressure  produced  by  closure  is  assumed  to  be  given 
by  the  relation:  p  =  ipcvj,  where  p  and  c  are  the  density  and  speed  of 
sound,  respectively,  and  vj  is  the  (relative)  impact  velocity  at  closure. 
This  result  was  first  derived  in  studies  of  the  water  hammer  effect  [  1 8 3 . 

Two  rectangularly-shaped  pressure  waves  are  produced  which  propagate  away 
from  the  closure  path  [1,15]-  However,  the  magnitude  and  directions  of 
these  closure  waves  are  also  influenced  by  the  time  at  which  closure  occurs. 
Time  of  closure  data  (also  obtained  from  the  second  stage  solution)  is 
plotted  in  Figure  5- 

Closure  is  predicted  to  first  occur  253-4  milliseconds  after  the  blast 
at  a  horizontal  range  of  725  ft.  Subsequent  closures  along  the  closure  path 
can  be  visualized  as  two  points  of  closure  moving  in  opposite  directions 
along  the  closure  path  away  from  the  point  of  first  closure.  The  horizontal 
speeds  of  these  points  are  given  by  the  dashed  curve  in  Figure  5-  Close  to 
the  initial  closure  point  the  speeds  of  both  points  of  closure  decrease 
rapidly  from  infinity.  The  horizontal  speed  of  the  outwards  moving  closure 
point  gradually  approaches  the  speed  of  sound  (5000  ft. /sec.)  until  it 
reaches  the  maximum  cavitation  range  567-5  milliseconds  after  the  blast,  but 
its  speed  is  always  supersonic.  The  inwards  moving  closure  point  drops 
rapidly  to  subsonic  speeds.  In  this  case  it  slows  to  less  than  700  ft.  per 
sec.  at  horizontal  ranges  near  150  ft.  from  the  blast  before  it  accelerates 
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to  an  infinite  speed  again  above  the  blast.  Final  closure  of  the  cavitation 
region  is  predicted  to  occur  directly  above  and  681  milliseconds  after 
the  blast.  (Effects  of  explosion  bubble  motions  are  also  ignored  in  this 
multi-stage  analysis.  For  an  explosion  of  this  magnitude,  bubble  pulse 
effects  are  expected  to  become  significant  about  one  second  after  the  blast[l6].) 

The  effect  of  rapid  changes  in  speed  of  the  closure  points  is  to 
produce  pressure  waves  with  curved  wavefronts.  The  effect  of  finite,  but 
supersonic,  closure  point  speeds  is  to  produce  pressure  waves  which  prop¬ 
agate  obliquely  to  the  closure  path.  Both  of  these  effects  can  be 
accounted  for  by  the  graphical  method  for  solving  the  third  stage  of  the 
cavitation  problem  which  is  illustrated  in  Figures  6  and  7. 

The  solid  line  (A)  in  Figure  6  shows  a  portion  of  the  closure  path 
near  the  region  of  first  closure,  but  is  now  plotted  with  equal  horizontal 
and  vertical  scales.  The  lower  dashed  line  (B)  in  the  figure  shows  a  portion 
of  the  closure  time  plot.  (The  rectangular  boxes  mark  computed  values  of 
closure  times, Tc,  which  are  also  tabulated  in  Appendix  A.) 

An  initial  wavefront(C)  can  be  constructed  from  the  closure  path  and 
time  curves  if  one  assumes  a  ray  theory  approximation  to  the  pressure  wave 
propagation  problem  [19].  To  construct  this  wavefront,  an  incremental  time,  T j , 
after  inita)  closure  is  first  selected  (in  this  case  T;  =  259*5  milliseconds). 

The  location  of  the  initial  wavefront  at  each  data  point  along  the  closure 
path  curve,  is  then  determined  by  drawing  a  vector  perpendicular  to  the  closure 
path  with  a  length,  L  -  c(Tc  -  Tj)  for  each  horizontal  range  where  closure 
occurs  before  Tj.  (A  second  wave,  not  shown  in  the  figure,  also  propagates 
downwards  from  the  closure  path.)  By  fitting  an  arc  to  the  initial  wavefront, 
its  center  of  curvature  (0)  is  graphically  determined.  The  ray  theory 

m 

approximation  is  generally  valid  if  the  resultant  radius  of  curvature  vector 
is  at  least  lOx's  longer  than  the  largest  vector  used  to  construct  the 


initial  wavefront. 


Subsequent  propagation  of  the  closure  pulse  pressure  wave  Is  deter¬ 
mined  by  constructing  a  series  of  arcs  of  concentric  circles  as  shown  in 
Figure  7*  The  arcs  are  constructed  about  the  initial  center  of  curvature 
until  they  intersect  with  the  water  surface.  Then  they  are  drawn  about  the 
surface-reflected  image  of  that  point.  Similar  wavefront  diagrams  can  be 
made  along  all  portions  of  the  closure  path  where  the  closure  point  speed 
is  supersonic  by  repeating  the  method  described  above  with  different 
initial  wavefronts.  An  example  of  a  region  farther  from  the  point  of  first 
closure  is  also  given  in  Figure  7- 

The  duration  and  magnitude  of  the  closure  pulse  at  an  arbitrary  point 
are  found  by  plotting  rays  perpendicular  to  the  constructed  wavefronts, 
from  that  point  back  to  the  closure  path.  For  example,  the  duration  of  the 
closure  pulse  at  a  horizontal  range  of  800  feet  and  a  depth  of  21  feet  is 
determined  by  the  difference  in  length  between  the  two  ray  vectors  (A)  and 
(B)  drawn  in  Figure  7.  The  magnitude  of  the  closure  pulse  is  given  by  the 
formula:  P]  =  pc  -  (r  j/ ( r  j  +  ra)),  where  pc  =  jpcv;,  ra  is  the  length  of 
vector  (A)  and  pj  is  the  peak  pressure  at  the  point  of  interest. 

In  order  to  permit  qualitative  comparisons  between  the  two  approaches, 
the  complete  pressure  history  predicted  by  the  multi-stage  approach  at  a 
horizontal  range  and  depth  of  800  ft.  and  21  ft.  respectively,  for  an 
explosion  of  A0000  lbs.  of  HBX  at  a  depth  of  200  ft.  is  shown  in  Figure  8. 


(13) 


5-  CONCLUSIONS 

Predictions  of  the  pressure  histories  of  bulk  cavitation  closure 
pulses  have  been  obtained  by  two  different  approaches.  A  hydrodynamic, 
finite-difference  method  directly  produces  pressure  histories  in  forms 
which  can  easily  be  compared  with  experimental  results.  However,  computa¬ 
tional  requirements  of  this  method  are  too  great  to  permit  its  convenient 
use  in  solving  the  full  2-0  problem.  (The  1-D  finite  difference  code 
used  in  this  paper  runs  in  about  the  same  time  --  1  minute  --  on  a  DEC 
10  computer  as  the  entire  2-D  multistage  code.)  A  2-D  multi-stage  method 
is  capable  of  characterizing  many  bulk  cavitation  features  over  the  entire 
region  of  interest,  but  the  approximate  formulations  on  which  the  solutions 
are  based  could  be  neglecting  significant  physical  mechanisms  or  could 
restrict  its  applicability. 

For  example,  the  need  to  extend  the  multi-stage  method  to  account 
for  retardation  of  the  bottom  boundary  motion  (lower  layer  accretion)  was 
discovered  by  comparing  results  of  the  two  different  approaches  [173- 
A  case  which  reveals  the  limited  applicability  of  the  multistage  approach 
is  illustrated  in  Figure  3  (which  was  obtained  by  the  finite  difference 
approach).  In  that  example,  a  bottom  reflected  shock  wave,  which  could  be 
described  by  a  first-stage  solution,  entered  the  cavitation  region  after 
the  second  stage  of  the  problem  had  commenced.  It  is  not  clear  to  this 
author  how  the  program  could  be  modified  to  include  such  cases  where  two 
stages  of  the  problem  overlap. 

It  should  be  noted  that  the  formula  used  in  the  graphical  method 
described  above  to  predict  the  magnitude  of  the  closure  pulse  is  not  con¬ 
sistent  with  previous  suggestions  by  other  investigators  [15].  This  new 

* 

result,  which  follows  directly  from  solutions  of  the  acoustic  wave  equation 
with  appropriate  boundary  conditions,  predicts  that  obliquely  propagating 


OM 

plane  closure  waves  will  have  significantly  higher  pressures  than  previously 
suggested. 

A  major  weakness  of  present  theoretical  methods  for  predicting  closure 
pulse  characteristics  is  the  inability  to  solve  cases  where  closure  occurs 
at  subsonic  speeds.  A  region  of  subsonic  closure  nearly  always  exists 
between  the  point  of  first  closure  and  the  blast.  In  the  example  used  in 
this  paper  the  subsonic  closure  region  extends  inwards  from  a  horizontal 
range  of  550  ft.,  which  is  175  ft  from  the  range  of  first  closure  --  less 
than  the  length  of  a  typical  Naval  vessel.  Anomalously  high  closure 
pulse  pressures  have  been  experimentally  observed  which  are  suspected  of 
being  effects  of  still  unidentified  subsonic  closure  mechanisms  [20]. 

Another  major  weakness  is  the  neglect  of  effects  due  to  pulsations  of  the 
vapor  bubble  produced  by  the  explosion. 

Further  improvement  of  numerical  bulk  cavitation  closure  pulse  pre¬ 
diction  methods  could  be  achieved  if  the  basic  mechanisms  of  bulk  cavitation 
were  better  understood.  However,  the  theoretical  2-D,  finite-difference 
approach  discussed  above  shares  common  disadvantages  with  experimental 
investigations  conducted  at  sea:  they  are  too  expensive  to  conveniently 
repeat  often  enough  to  explore  parameters  of  the  problem.  These  para¬ 
metric  studies  are  needed  to  first  identify  and  then  quantify  new  cavitation 
mechanisms,  such  as  those  which  might  be  controlling  subsonic  closure 
effects.  Laboratory  facilities  for  generating  easily  repeatable,  small 
scale  bulk  cavitation  effects,  therefore,  appear  to  provide  a  critically 
needed  tool  for  studying  bulk  cavitation  mechanisms. 
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FIGURE  (I):  RESPONSE  OF  A  PRESSURE  GAGE  TO  AN  UNDERWATER 
EXPLOSION  SHOWING  A  CLOSURE  PULSE  PRODUCED  BY 
BULK  CAVITATION. 
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SURFACE  VELOCITY  AND  PRESSURE  HISTORIES  ABOVE  AN  UNDERWATER 
EXPLOSION  SH0WIN6  EFFECTS  OF  A  REFLECTINS  BOTTOM. 
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FI6URE  (5):  TIME  OP  CLOSURE  (A)  AND  HORIZONTAL  SPEED  OF  CLOSURE  (B)  OF 
THE  CAVITATION  REGION  PLOTTED  AS  A  FUNCTION  OF  RANGE 
FROM  AN  UNDERWATER  EXPLOSION.  RESULTS  PREDICTED  BY 
MULTI-STAGE  METHOD. 
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FIGURE  (6):  CLOSURE  PATH  (A)  AND  TIME  CB)  DATA  USED  TO  GRAPHICALLY 
DETERMINE  AN  INITIAL  WAVEFRONT  CC)  OF  THE  CLOSURE  WAVE 
AND  ITS  CENTER  OF  CURVATURE  CD). 
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FIGURE  (7):  WAVEFRONT  AND  RAY  CONSTRUCTIONS  USED  TO  GRAPHICALLY 
DETERMINE  THE  CLOSURE  PULSE  FORM  AT  AN  ARBITRARY 
POINT  ABOVE  THE  CLOSURE  PATH. 
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FIGURE  (8):  PRESSURE  HISTORY  AT  A  POINT  21  FEET  DEEP  AND  880 
FEET  FROM  AN  UNDERWATER  EXPLOSION  AS  PREDICTED  BY 
THE  MULTISTAGE  METHOD. 
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